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The density of states of graphene has Van Hove singularities that can be reached by chemical 
p^ I doping and have already been explored in photoemission experiments. We show that in the presence 

^^ , of Coulomb interactions the system at the Van Hove filling is likely to undergo a Pomeranchuk 

^v^ ■ instability breaking the lattice point group symmetry. In the presence of an on-site Hubbard 

interaction the system is also unstable towards ferromagnetism. We explore the competition of 
the two instabilities and build the phase diagram. We also suggest that, for doping levels where 
. ■ the trigonal warping is noticeable, the Fermi liquid state in graphene can be stable up to zero 

^_^ ' temperature avoiding the Kohn-Luttinger mechanism and providing an example of two dimensional 

^^^ \ Fermi liquid at zero temperature. 

^ ^ ■ PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 

I. INTRODUCTION 

The recent synthesis of single or few layers of graphite [1,|2| has permitted to test the singular transport properties 
predicted in early theoretical studies^, [j,[p,Q and experiments [TJ. The discovery of a substantial field effect [8,] and 
the expectations of ferromagnetic behavior Q has risen great expectations to use graphene as a reasonable replacement 
of nanotubes in electronic applications. 

^O ■ The Fermi level of graphene can be tuned over a wide energy range by chemical doping [10, lU , or by gating P, M\ • 

C I The full dispersion relation of graphene (Fig. 2) shows several regions of special interest. The low doping region near 
^ . half filling has been the object of much attention due to the behavior of quasiparticles as massless Dirac electrons. 
The Fermi surface at low filling (cusps of the dispersion relation) begins being circular. For increasing values of the 
electron density the trigonal distortion appears and, finally, several Van Hove singularities (VHS) develop at energies 

Cn ■ the order of the hopping parameter E ~ 2.7ey. The possibility of finding interesting physics around these densities 
^ , have been put forward recently in an experimental paper reporting on angle resolved photoemission (ARPES) results 
[l^. Chemical doping of graphene with up to the VHS has been demonstrated in [10]. The VHS are of particular 
interest in the structure of graphene nanoribbons [12l, 13] and nanotubes where observation of Van Hove singularities 

l/^ ■ using scanning tunneling microscopy have been reported [14]. It is known that around the Dirac point due to the 
vanishing density of states at the Fermi level, short range interactions such as an on-site Hubbard term U are irrelevant 
in the renormalization group sense [15, 16] and should not give rise to instabilities at low energies or temperatures. 
The possibility to stabilize a ferromagnetic [13] or superconducting [l^, ^^, [20] phase near half filling is very unlikely. 

(^ , Alternatively, at densities around the VHS the physics is dominated by the high density of states and should resemble 
the one discussed before in the framework of the high- Tc superconductors [2l| . This is the physical situation to look 
for ferromagnetic or superconductivity [22[ in graphene -or graphite-. In ref. [23| it was argued that the most likely 
instabilities of a system whose Fermi surface has VHS and no special nesting features are p-wave superconductivity 
^ I and ferromagnetism. A very interesting competing instability that may occur when the Fermi surface approaches 
singular points is a redistribution of the electronic density that induces a deformation of the Fermi surface breaking 
the symmetry of the underlying lattice. This phenomenon called Pomeranchuk instability ^] has been found in the 
squared lattice at the VHS filling [25|, [20, [23, [28| and has played an important role in the physics of the cuprate 
superconductors [29, 30] and in general in layered materials [31]. It has also been discussed in a more general context 
in [32] where it was argued that opening of an anisotropic gap at the Fermi surface is an alternative to cure the 
infrared singularities giving rise to Pomeranchuk instabilities. More recently the Pomeranchuk instability has been 
studied in relation with the possible quantum critical points in strongly correlated systems [33]. Due to the special 
symmetry of the Fermi surface of the honeycomb lattice, the pattern of symmetry breaking phases can be richer than 
these of the square lattice. 

In this paper we study the Van Hove filling in graphene modelled with a single band Hubbard model with on-site U 
and exchange Coulomb interactions V. We perform a mean field calculation and show that Pomeranchuk instability 
occurs very easily in the system in the presence of the exchange interaction V. Adding the on-site interaction U allows 
for ferromagnetic ground states. We study the coexistence or competition of the two and build a phase diagram in the 
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FIG. 1: (Color online) The honeycomb lattice of graphene. 



(U, V) space of parameters. In section |TT] we construct the exchange V and on-site U interactions in the honeycomb 
lattice and describe the method of calculation and the results obtained. In the last section we discuss some points 
related with the physics of graphene at high doping levels. In the context of the original Kohn-Luttinger instabilities 
[34] that establishes that all Fermi liquids will be unstable at low enough temperatures, we suggest that graphene at 
intermediate fillings can provide an example of stable metallic system at zero temperature ^, i36|, [33] . Finally we set 
the lines for future developments. 

II. THE MODEL 

The deformation of the Fermi surface by the electronic interactions is a classical problem in condensed matter dating 
back to founding work of Luttinger [38] . We model the system by a single band Hubbard model in the honeycomb 
lattice with an on-site interaction U and a nearest neighbor Coulomb interaction V, and perform a self-consistent 
calculation along the lines of refs. [261,1391. 

The Hubbard hamiltonian is 

a = t^Y^cf^c\a + t/^niimi +V ^ ni^nj^/ . (1) 

i;cr i <ij;crcr'> 

The hopmng parameter t of graphene is t ~ 2.1 eV . The next-to-nearest neighbor hopping t' is estimated to be much 
smaller [40| and it will be ignored. Due to the special topology of the honeycomb lattice with two atoms per unit cell 
the Hamiltonian is a 2 x 2 matrix. We will work out in some detail the different terms. The operators a, h associated 
to the two triangular sublattices A, B are 

where N a = ^b = A/" is the number of lattice cells and ^i are the three vectors connecting a point of sublattice A 
with its three neighbors in sublattice B (see Fig. [1]). 
In term of these operators the free Hamiltonian reads 



ka 



where we have included a chemical potential /i. 
In matrix form the non-interacting Hamiltonian is 



t(/)*(k) /i 
where 
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FIG. 2: (Color online) Left: Lower band of the dispersion relation of graphene. The black line marks the Fermi surface for the 
Van Hove filling. Right: Evolution of the Fermi surface of graphene with doping. At half filling the Fermi surface consists on 
the six vertices of the hexagonal Brillouin zone (black points). Only two are independent, the rest being equivalent by lattice 
vectors. Emptying the half filled band the Fermi surface develops circular hole pockets that undergo the trigonal distortion (in 
red) and become the yellow triangles at the Van Hove fillig jjl — t. 



The Hamiltonian (|4]) gives rise to the dispersion relation 



V3 



V3 



e^{k) =11^ ty |^(^)p = /i ± ty 1 + 4cos^ — akx + 4 cos — akx cos -aky , 



whose lower band is shown in the left hand side of Fig. [2l 



III. COULOMB INTERACTION V AND POMERANCHUK INSTABILITY 



(6) 



We will begin by studying the influence of the interaction V omitting the spin of the electrons. The on-site U will 
be added to study the competition of the Pomeranchuk instability with ferromagnetism. The interaction V takes 
place between nearest neighbors belonging to opposite sublattices. It reads: 



Hv = V^ [ai^ak+q6i^/6k/-q0*(q) + 6i^6k+qai^,ak/-q^(q)]. 

kk'q 

The mean field Hamiltonian becomes 

Hmf = ^[^(k)a+6k + ^*(k)6+ak] + /i ^{ci^cti^ + ^k ^k), 



(7) 



(8) 



where 



V 



<^(k)=#(k)-^T/*(q)<&++qak+q> , V{ci) = -J2r{cD 



(9) 



and A'' is the number of lattice cells. In matrix form and including the chemical potential the mean field Hamiltonian 
reads: 



H. 



M t<P{k) - Eq V* (q) < b+^^ak+g > 



The mean field Hamiltonian (p!Q|) gives rise to the dispersion relation 



(10) 



(11) 




FIG. 3: Pomeranchuk instability suffered by an initial Fermi surface at the Van Hove singularity (upper part) when an exchange 
interaction V=lt is added (lower part). 



where F(y, k) is 



+v^ Eqq' r (q)</'(q') < 4+gh+q >< ht+^,ak+q' > . 



(12) 



We look for a self-consistent solution of eq. (pT]) imposing the Luttinger theorem i. e. that the area enclosed by 
the interacting Fermi line is the same as the one chosen as the initial condition. We begin with an initial (free) Fermi 
surface, as an input, add the interaction and let it evolve until a self-consistent solution is found with a given "final" 
interacting Fermi surface. The important constraint is that the total number of particles should remain fixed. This 
procedure has been used in the same context in [261, [39] . 

The Fermi surface of graphene around the Van Hove singularity is shown in Fig. [2l To better visualize the 
Pomeranchuk instability we keep the image of the neighboring Brillouin zones. Fig. [3] shows the initial Fermi surface 
sitting at the VHS (upper part) and the spontaneous deformation obtained self-consistently from eq. (pT]) for a value 
of the exchange interaction V = t. 

IV. FERROMAGNETISM AND COMPETITION OF FERROMAGNETISM AND POMERANCHUK 

INSTABILITIES 

In order to study a possible ferromagnetic instability we will add to the free Hamiltonian an on-site Hubbard term 

Hu = U^lat^ai^al^aii -^ bt^bi^bl^bii], (13) 



A mean field ferromagnetic state will be characterized by 

with 

I n m _ I _ n m 

what produces the mean field Hamiltonian 



(14) 
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FIG. 4: (Color online) Left: Initial polarization of the Fermi surface with the spin up electrons (upper part) less populated 
than the spin down (lower part): ni — n^ = 0.14, /i| = — 0.9t, /x^ = — l.lt. Center: Final state obtained for the values of the 
parameters U=lt and V=2t. A Pomeranchuk instability is clearly visible in both spin bands. Right: Final state with enhanced 
ferromagnetic polarization obtained for the values of the parameters U=2t and V=lt. 



Adding the kinetic term we get the Hamiltonian 
whose dispersion relation 



/i + f(n + crm) t(j){k) 

t(^*(k) /i + f(n + crm) 



^'/ca =Ma ±^|^(k)| , 11^= 11^ —{n^crm), 



(17) 



(18) 



is solved self-consistently to get the ferromagnetic ground state. 

The competition of the Pomeranchuk instability found in the previous section with a possible ferromagnetic in- 
stability is studied with the same procedure using the full mean field Hamiltonian Hq + H^^^ + ^mf ^^^ solving 
self-consistently the equations 



^T(k) = M + f/n^ - ./t2|^(k)|2 + j^(y^k)^, 



(19) 



Ei(k) =^,+ Un^- ^t2|^(k)|2 + F(F,k)x, 
where F(V, k)cr is given in eq. (12) and Ucr = l/2(n + am). 



(20) 



V. SUMMARY OF THE RESULTS 



The competition of the Pomeranchuk instability with ferromagnetism is exemplified in Fig. [H The figure in the 
left side represents the initial free Fermi surface where the spin up (upper side) and down (lower side) electrons have 
different populations around the Van Hove filling. In particular in the example given, the chemical potential of the 
spin up (down) electrons is set slightly below (above) the VHS: ni — n^ = 0.14, ji^ = — 0.9t, /i| = —l.lt. The image 
in the center represents the renormalized Fermi surface for the two spin polarizations when the exchange interaction 
V is bigger than IJ: V = 2U = 2t. We can see that the ferromagnetism has disappeared: the final Fermi surface is 
the same for the two spin polarization and presents a Pomeranchuk deformation. The opposite case is shown in the 
figure at the right: with the same initial configuration the final state for the values of the interactions U = 2V = 2t 
has an enhanced ferromagnetism with no signal of deformation. 
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FIG. 5: Phase diagram showing the competition of Pomeranchuk and ferromagnetic instabihties as a function of the interactions 
U and V measured in units of the hopping parameter t. 



The phase diagram of the dominant instabihty as a function of the interaction strength U and V (measured in units 
of the hopping parameter t) is shown in Fig. [5l In all cases the initial Fermi surface consists of a slightly polarized 
state around the Van Hove singularity as the one shown in the left hand side of Fig. [H The symbols in Fig. [5] denote 
calculated points with the following meaning: the circles appearing for low values of U denote an unpolarized final 
state with Pomeranchuk deformation as the one shown in the center of Fig. [H Crosses appearing for large values of 
U represent a final state where the polarization is bigger than the initial one and the Fermi surface has the original 
symmetry as the one in the right hand side of Fig. [H The points denoted by correspond to values where both 
instabilities coexist and the final Fermi surface is deformed and spin polarized. Finally the crosses at low values of U 
with V=0 represent final states which are still polarized but where the spin polarization is smaller than the initial one. 
This region is denoted by FM | in the figure phase diagram. PI (FM) denotes the region where the Pomeranchuk 
(ferromagnetic) instability dominated respectively. 

The ferromagnetic evolution of the system for V=0 is as follows: Below a critical value of 1.5 < /7c < 2 a free 
slightly polarized system at fillings near (but not at) the VHS evolves towards an unpolarized one in agreement with 
[4l|, [42|- The Coulomb exchange V induces a deformation of the Fermi surface already at very small values Vc ~ OAt 
when the free initial state is near the Van Hove singularity. Increasing U enhances the magnetic polarization of the 
interacting system until the critical value U ~ 1.75t is reached where the final state corresponds to a ferromagnetic 
system with un-deformed Fermi surface. Unlike what happens in the square lattice, the critical value of U above 
which ferromagnetism prevails does not depend on V (vertical line in Fig. [5]). In the blank region in the upper part 
of the figure corresponding to high values of V around the critical U we have not been able to reach a self-consistent 
solution. 

We note that at V=0, the critical U for ferromagnetism is zero at the Van Hove filling and changes very rapidly 
around it in a rigid band model [4l|, [42[ . This behavior is due to the divergent density of states at the VHS that 
would be smeared by temperature effects or disorder in real samples. Based on the previous studies if the square 
lattice we are confident that the results obtained in the present work are robust and the phase diagram shown in fig. 
[5] will remain qualitatively when more detailed calculations are done. 

The realization of one or another phase in real graphene samples depends on the values that the effective Coulomb 
interaction parametrized by [/, V has at the VHS filling. Besides the ARPES experiments [10] which do not comment 
on the strength of the interaction, we are not aware of other experiments at these doping levels. Conservative estimates 
can be obtained from the graphite intercalated compounds [43] on the basis of the similarities found in [W\ between 
the Fermi surfaces of the two. Values of [/ ~ 2.5t, V ^ t can be very reasonable and lie in the range discussed in the 
present work. 



VI. DISCUSSION AND OPEN PROBLEMS 

In this work we have examined some physical issues to be expected in graphene at high doping. We have seen that 
at the Van Hove fihing in the presence of an exchange Coulomb interaction the Fermi surface is softer and becomes 
very easily deformed as compared with a similar analysis of the square lattice. 

A natural extension of the work presented in this article is to study the competition of the instabilities studied 
in this work with superconductivity along the lines of ref. [44] . A complete renormalization group analysis of the 
competing low energy instabilities of the Van Hove filling as a function of the doping and the couplings U and V is also 
to be done. Due to the special geometry of the Fermi surface in the honeycomb lattice the symmetry breaking phases 
can present a richer variety than those of the square lattice. A fairly complete analysis of the short range interactions 
around the Dirac point was done in the early paper [18]. Due to the existence of two Fermi points the classification 
of the possible low-energy couplings is similar to the g-ology of the one dimensional models. The RG classification 
of the low-energy couplings in the Van Hove filling is richer since there are three independent VH points to consider 
and work in this direction is in progress. An analysis of the physics of the Fermi surface around the trigonal warping 
can lead also to very interesting results in the light of the evolution of the Fermi surface anisotropics done in |45|, l4y] . 

On the view of the analysis of this problem done in the square lattice we can expect that RG calculation will enrich 
the phase diagram but the two phases discussed here will stay. Ferromagnetism is a likely possibility due to the Stoner 
criterium although it will compete with antiferromagnetism at high values of U. A very important parameter to this 
problem is the next to nearest neighbor hopping t' that suppresses the nesting of the bare Fermi surface and affects 
the shape of the phase transition lines [4l|, [42|. These references focuss on the magnetism of the Hubbard model 
(U) on the honeycomb lattice at finite dopings and found that non-homogeneous (spiral) phases are the most stable 
configurations around the Van Hove filling. It would be interesting to analyze the competition and stability of these 
non-homogeneous ferromagnetic configurations in the presence of a V interaction. 

Of particular interest is the possibility of coexistence or competition of the Pomeranchuk instability with possible 
superconducting instabilities [22]. A very interesting suggestion has been made in studies of the square lattice that 
superconductivity changes the nature of the Pomeranchuk transition going from first to second order [4J] . This feature 
will probably be maintained in the honeycomb lattice and we are actually exploring this problem. 

As it is known, the Kohn-Luttinger instability in its original context [3J] suggests that no Fermi liquid will be 
stable at sufficiently low temperatures. In the case of a 3D electron system with isotropic Fermi surface there is 
an enhanced scattering at momentum transfer 2/cf which translates into a modulation of the effective interaction 
potential with oscillating behavior similar to the Friedel oscillations. This makes possible the existence of attractive 
channels, labelled by the angular momentum quantum number. The issue of the stability of Fermi liquids was studied 
rigorously in two dimensions in a set of papers [35, ^, 37] with the conclusion that a Fermi liquid could exist in 
two space dimensions at zero temperature provided that the Fermi surface of the system obeys some " asymmetry" 
conditions. In particular the non-interacting Fermi surface had to lack inversion symmetry 

e(-k) ^ e(k) 

at all points k and be otherwise quite regular (the Van Hove singularities can not be present). The dispersion relation 
of graphene obeys such a condition for fillings where the trigonal warping is already noticeable and below the Van Hove 
filling provided that disorder and interactions do not mix the Van Hove points. This is a very common assumption 
in the graphene physics around the Dirac points where the inversion symmetry in k-space acts as a time reversal 
symmetry [47]. We find interesting to note that the -quite restrictive-conditions of refs. [35|, [3y, [33 can be fulfilled 
in a real -and very popular-system. 
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